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Abstract 

', Group classification of systems of two coupled nonlinear reaction-diffusion 

O I equation with a diagonal diffusion matrix is carried out. Symmetries of diffusion 

^ I systems with singular diffusion matrix and additional first order derivative terms 

^ ■ are described. 

^ ; 
O ■ 

; 1 Introduction 
in '■ 

>■ ■ Coupled systems of nonlinear reaction-diffusion equations form the basis of many 

^ . models of mathematical biology. These systems are widely used in mathemat- 

O ! ical physics, chemistry and also in social sciences and many other fields. Such 

^ I reach spectrum of applications stimulates numerous thorough investigations of 

I fundamentals of these equations theory. 

O ■ In the present paper we continue group classification of systems of reaction- 

^ . diffusion equations with general diffusion matrix 

J^' Ut-Am{A^^u + A^h)=f\u,v), 

t^: vt-AUA''u + A''v) = f{u,v) 



where u and v are function of t,Xi,X2, ■ ■ ■ ,Xm, A^^ , A^"^ , A"^^ and A"^"^ are real 
constants and is the Laplace operator in K^. 
^ ' Up to linear transformations of functions u, v and p it is sufficient 

to restrict ourselves to such diffusion matrices (i.e., matrices whose elements 
are A^^,--- ,v4^^) which are diagonal, triangular, or are sums of the unit and 
antisymmetric matrices. In the last case ([T]) can be reduced to a single equation 
for a complex function ( generalized complex Ginzburg-Landau (CGL) equation) 
whose group classification was carried out in paper pLj. 

In the present paper we classify equations ([I]) with a diagonal diffusion matrix. 
Without loss of generality such equations can be written as 



vt - aA„,v = P{u,v) 
where a is a constant. The related diffusion matrix is A 



(2) 
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Just equations of type are the most popular models of reaction-diffusion 
systems first introduced by Turing in 1952 [2]. It is practically impossible to 
enumerate all fields of applications of such equations. We restrict ourselves to 
few examples only. 

• The complex Ginzburg-Landau (CGL) equation 

Wt- {l + i(3)A2W = W - {l + ia)\W\'^W (3) 

can be presented as a system ([T]) where u and v are the real and imaginary 
part of W. In particular case (3 = this system takes the form (l2l). 

• The primitive predator-prey system which can be defined by [12] 

Ut - DUxx = -uv, .^s 
Vt — \DVxx = uv 

also appears as an particular subject of our analysis. 

• The X — uj react ion- diffusion system [13] 

Ut = A2U + X{R)u - uj{R)v, . s 

Vt = A2V + uj{R)u + X{R)v, 

where R^ = u"^ + v^, is widely used in studies of reaction-diffusion models, 
in particular, to describe spiral waves phenomena [H] . 

Symmetries of equations were studied in paper ^3]. We shall add the 
results [15] in the following. 

• The Jackiw-Teitelboim model of two-dimension gravity with the non-relativistic 
gauge [To] appears as a particular ((l-l-l)-dimensional) case of the following 
system: 

Ut — AmU = 2ku — 2u^v = 0, , . 

Vt + AmV = 2uv^ -2kv = 0. ^ ^ 



Symmetries of equations ([6]) for m = 1 were investigated in paper pi]. In 
the following we complete the results obtained in [TT] . 

Apparently the first attempt of group classification of equations ([2]) was made 
by Danilov [1]. But the results present in [3] are rather incomplete. 

Group classification of equations with general non-degenerated diffusion 
matrix was announced in [3] and presented in [U]. However, the equivalence 
relations where not used systematically there to simplify the equations which 
resulted in rather cumbersome form of the classification results. Moreover due to 
typographical errors the tables with classification results present in [6] are poorly 
readable (see [T] for additional comments). 
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Group classification of systems of heat equations 

ut - {f{u,v)u^)^ = 0, 
vt - {g{u,v)v^)^ = 

has been performed in paper [7]. For some classes of functions f{u, v) and g{u, v) 
equations ([7j) are equivalent to ([2]) and the related classification results can be 
compared with ones obtained in the previous section. We will do it in Section 8. 

Symmetries of systems of reaction-diffusion equations with a diagonal diffu- 
sion matrix (i.e., of systems (E])) where studied in papers [8], [9]. We will show 
in the following that the classification results obtained in [8], [9] are incomplete 
and include many equivalent cases treated as non-equivalent ones. 

The problem of group classification of equations is still relevant and we 
will present its solution here. In addition, we classify equations with non- 
invertible diffusion matrix (i.e., equations (E]) when parameter a is equal to zero) 
and also the following equations with first order derivative terms: 

where Ux^ = Pi_i are arbitrary constants and summation from 1 to m is 
imposed over the repeated index Moreover, without loss of generality one can 
set 

Pl=P2 = ■■■ = Pm~l =0, Pm=P, P = 0, 1- (9) 

In the case p = equation ([H]) reduces to ([2]) with a = 0. We notice that 
this equation is used in such popular models of mathematical biology as the 
FitzHung-Naguno [16] and Rinzel-Keller [TTj ones. 



(7) 



(8) 



2 Equivalence transformations 

The problem of group classification of equations ([2]) and ([8]) will be solved up to 
equivalence transformations. Clear definition of these transformations is one of 
the main points of any classification procedure. 
We say the equations 

Vt - aA^v = P{u,v) 

be equivalent to ([2]) if there exist an invertible transformation u —>■ u = G{t, x, u, v), 
V ^ V = $(t, x,u,v), t ^ t = T(t, x,u,v), X ^ X = X(t, X, u, v) and ^ = 
F^i^UjtjX, p, P) which connects ([2]) with ( fTOl) . In other words the equivalence 
transformations should keep the general form of equation ([21) but can change the 
concrete realization of non-linear terms and 
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The group of equivalence transformations for equation can be found using 
the classical Lie approach and treating and p as supplemental dependent 
variables. In addition to the obvious symmetry transformations 

t^t' = t + a, x^—>x'^ = R^uXu + h^ (11) 

where a, and R^j^u are arbitrary parameters satisfying R^^R^x = S^\, this group 
includes the following transformations 

t ^ \~H, Xb ^ \-^Xb ^ ^ 

and 

Ub ^ K^'u,, /" a-'k^'^r, .^3. 
t — > a^^t, Xb ^ Xb, a ^ 

where indices b, c take values 1 and 2, K^'^ and K^"^ are elements of invertible 
constant matrices K and K respectively, moreover, K commutes with A, and K 
satisfy the condition 

kA{a)K-^ = aA{l/a); 

A 7^ and ba are arbitrary constants, and we use the temporary notations u = 

Ui,V = U2- 

If parameter a is equal to 1 then K is an arbitrary 2x2 invertible matrix, and 
equivalence transformations ( |T3l) is trivial. If a 7^ 1 then K and K are arbitrary 
non-degenerated diagonal and anti-diagonal matrices respectively . 

Transformation (1131) reduce to the change a — >■ 1/a in the related matrix 
A, i.e., to scaling the free parameter a. Thus without loss of generality we can 
restrict ourselves to the following values of a: 

1. a = 0, 

2. - 1 < a < 0, < a < 1, (14) 

3. a = 1. 



It is possible to show that for a 7^ there is no more extended equivalence 
relations valid for arbitrary nonlinearities and p. If a = there exist a pow- 
erful equivalence relation m — m, f — '^{v) with an arbitrary function ^p{v). 
However for some particular functions p and p the invariance group can be 
more extended. In addition to transformations (fT2|) and (fT3|) it includes symme- 
try transformations which does not change the form of equation Moreover, 
for some classes of functions P,P equation ([2D admits additional equivalence 
transformations (AET) which belong neither to symmetry transformations nor 
to transformations of kind (fT2l) . ( fT3l) . 
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In spite of the fact that we search for AET after description of symmetries 
of equations and specification of functions f^, p, for convenience we present 
the hst of the additional equivalence transformations in the following formulae: 

1. M — s> e^^u, V — i> e''*f , 

2. u ^ u + ut, t> — > f , 

3. u u, V V + pt, 

4. u ^ u + fipt, V —>■ e^P^v, 

5. u e^^u, V ^ V — upt, 

6. M — > M, i; — > t> + ptu, 

7. u ^ e^'^^u, V V + ut"^, 

8. u ^ u + ut"^, V f e^'^*, 

9. u ^ u, V ^ V — 2ptu + pSt'^, 

10. u e^P^u, V e^^* {v + utu + pt'^u) , 

11. u ^ u + ript, V V — pt, 

12. u e'**M, V e'^^{v — uKtu), 

13. u^u + 2pt, V + 2ptu + 2pH'^, (15) 

14. u e'^'u, V e^^v, 

15. u e'^'^* {u cos{ujat) + v sin(co'crt)) , 

V e^"^* {v cos{uj(7t) — u sm{ujat)) , 

16. u e^"^* (m cos(crco't^) — v sin(crco't^)) , 

V e^"^* (-u cos(crcjt^) + u sm.{aujt^)) , 

17. u e^"^*' (m cos(2cjt) + V sin(2u;t)) , 

V e^"^*^ cos(2cut) - u sin(2co't)) , 

18. u e'^'^^u, V e'"^\v - autu), 

19. u e^^'^'u, V e^"^*' {v + 2utu), 

20. u e^'^*M, V — > e^'^*''t;, 

21. u^u + 3ujt, V + 3ujt^u + 3ujH^ + ptu + 3ujpt^, 

22. U ^ U + pXm, V ^ V. 

Here the Greek letters denote parameters whose values are either arbitrary 
or specified in the tables presented below. Equivalence transformations (ITSll are 
valid only for particular non-linearities which will be specified in the following. 



3 Symmetries and classifying equations 

We search for symmetries of equations ([2]) and (IHl) with respect to continuous 
groups of transformations using the infinitesimal approach. Applying the Lie 
algorithm or its specific formulation proposed in [6] one can find the determining 
equations for coordinates r/, vr^, vr^ of generator X of the symmetry group: 

X = r/9i + r5,.„-7rU-7r2a, (16) 
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and classifying equations for non-linearities and We will not reproduce 
the related routine calculations but present the general form of symmetry X for 
equation ([2]) with a 7^ 0, 1 [IJ : 

X = XK + a^G^" + cu^G'^ + jjD - C^udu - C^vd, - B^du - B^d^, , „^ 

where the Greek letters denote arbitrary constants, B^, B^ and C^, are func- 
tions of t, X and t respectively, and 

K = 2t{tdt + x^d^^) - ^{udu + \vd^) - tm{udu + vd^), 

D = tdt + lx^d,^, 

G^L = td^^-lxf,{udu + ^vdy), 

Gf, = e"^* (9^.^ - l-fXf,{udu + Ivd^)) ■ 

If a = then the related generator X again has the form (ITTll where however 
X = a'^ = uj^ = C'^ = and B^ is a function of t, x and u. 

For a = 1 the symmetry group generator (which we denote by X) is more 
extended and has the following form: 

X = X + C\d^ + C\du (19) 

where X is given in (fT7|l and C^, are functions of t. 

Equation ([2]) admits symmetry (fT9|) iff the following classifying equations for 
P and p are satisfied [1]: 

(A(m + 4)t + /i + |Ax2 + a^'x^ + 7e^*cu%) + C^) + C^p + C/m 

+ B] - AmB^ = {B^du + B^d., + C^udu + C^vd^ + C^ud^ + CSa„ 
+\mt{udu + t;<9^) + (|Ax2 + a^'x^, + 7eT'*cj%) (m9„ + Ivd^)) p, 



(A(m + 4)t + /i + i (iAx2 + a^^x^ + 76^*0;%) + C^) p + C^p + C^v 
+Cfu + 52 - AmB^ = (B^du + B^d^ + C^ud^ + C\d^ + C^ud^ + CS9„ 
+Amt(Ma„ + vd^) + (|Aa;2 + a^'x^, + 7e^*w%) (m(9„ + \vd^)) p. 

(20) 

In other words, to make group classification of systems (Ej) means to find all non- 
equivalent solutions of equations (l20l) and to specify the related symmetries (JTTIl 
[H]. We note that equations fpUl) can be decoupled equating terms multiplied by 
the same variables x^ or their powers. 

Consider now equation ([8]) and the related symmetry operator ( fT6l) . The 
determining equations for r], and vr*^ are easily obtained using the standard 
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Lie algorithm: 



r]tt = Vx^, = Vu = Vv = 0, 

= e = e = o, 

TT^ = TT^ = TT^ = nvr^ = 0- 

<,« + <,. = 0, ^''^ 

where subscripts denote derivatives w.r.t. the corresponding independent vari- 
able, i.e., Vt = ^, C = IE' etc. 

Integrating system (12T1) we obtain the general form of operator X: 

X = udt + pud,^ + ^"''d.x^ + pD- B'du - B^d, - Fudu - GvO,; (22) 
/i = 2(F - G) if p ^ (23) 

where B^, B^ are functions of {t,x), F and G are functions of t and summation 
over the indices a, v is assumed with cr, z/ = 1,2, ■■■,n — 1. 

The classifying equations for and p reduce to the following system 

(/i + F)f' + Ftu + (9t - A„)5i = [B'du + 5^9, + Fudu + /\ 

(24) 

(/i + 6-)/^ + Git; + Bl - pBl^ = {B'du + B^d, + Fud^ + Gvd,) f. (25) 

Solving (Ell), ( |25l) we shall specify both the coefficients of infinitesimal oper- 
ator fl22|) and the related non-linearities and 

It is obvious that the widest spectrum of symmetries corresponds to the case 
when the parameter a is equal to 1 since the corresponding generator X ffTIH) 
includes two additional terms C^ud^ and C^vdu- Quite the contrary, equations 
([2]) with a 7^ 1 and especially ([H]) admit relatively small variety of symmetries. 



4 Classification of symmetries 

Following [IJ we specify basic, main and extended symmetries for the analyzed 
systems of reaction-diffusion equations. 

Basic symmetries are nothing but generators of transformations (fTTi) forming 
the kernel of a symmetry group, i.e., 

Po = dt, P^ = d,^, J"" = x^d,^ - x,d,^. (26) 

Main symmetries form an important subclass of general symmetries (fT7j) and 
have the following form 

X = -pD + C^udu + C^vd^ + C^ud^ + C^vdu + B^d^ + B'^d^ (27) 
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(if a ^ 1 then C"^ = = 0). 

In accordance with the analysis present in pi| the complete description of 
general symmetries (fT71) can be obtained using the following steps: 



Find all main symmetries fl27p . i.e., solve equations (1201) for ^^^'^ 

{11 + + C^p + Clu + Cfv + Bl - A^B^ 

= {C^udu + C^vdy + C^ud^ + C^vdu + B^d^ + B^d^)f\ 

{fi + C^)p + C'p + C?v + C^u + 52 - aA^52 

= {C^udu + C^vd^ + C'^ud^ + C^vdu + B^du + B^d^)p. 



(28) 



Specify all cases when the main symmetries can be extended, i.e., at least 
one of the following systems be satisfied: 



ap = {audu + vd^)p 
p = [audu + vdv)f 



2_f...^ , „,;,^.2 (29) 



2 , , „,A^.2 (30) 



a{p + 7u) = (aM(9„ + vdy)p, 
p + ^v = {audu + vdy)f 

or if equation ( l29l) is satisfied together with the following conditions: 
(m + 4)/^ + /i^7i + 

= {{ij}^u + /i^^w + mu)du + (/i^^u + /i^^i' + mv)d^) p, 

(31) 

((Z/I^M + Z/12t;)9, + (z/21u + Z/^^i;)^,) 

were c = 1, 2, /i^^ and z/'^* are constants satisfying (a— l)/i^^ = (a— l)z/'^^ = 0. 

If relations fEUl) . flHUj) or flHTl) are valid then the system ([2]) admits symmetry 
or the conformal symmetry K — {tfi^^ + i'^^)udu — {tn^"^ + v^'^)vdu — 
(t/i^^ + iy'^^)udy — (t/i^^ + iy'^'^)vdy correspondingly. 

• When classifying equations ([HI) for p 7^ the second step in not needed 
since in accordance with (l22l) these equations admit only basic and main 
symmetries. 

In the following sections we find main and extended symmetries for the clas- 
sified equations. For clarity we start with group classification of systems (IHl) 
with p ^ which is more simple technically and present rather detailed calcula- 
tions. Then we consider equations ([2]) and present classification results without 
technical details. 
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5 Algebras of main symmetries for equation (IHl) 



To describe main symmetries we use the trick discussed in [T], i.e., make a priori 
classification of low dimension algebras of these symmetries. In accordance with 
(122|) any symmetry generator extending algebra (126!) has the following form 

X = fiD- B^du - B^d^ - Fudu + " ^) ^^v- (32) 

Let and X"^ be operators of the form ( |32l) then the commutator 
is also a symmetry whose general form is given by (|32l) . Thus operators ( |32l) form 
a Lie algebra which we denote as A. 

Let us specify algebras A which can appear in our classification procedure. 
First consider one-dimensional A , i.e., suppose that equation ([H]) admits the 
only symmetry of the form ( l32l) . Then any commutator of operator ( !26l) with 
( 132|) should be equal to a linear combination of operators ( l26l) and ( l32l) . Using 
this condition we come to the following possibilities only: 

X = X^ = HD- a^du - a^d^ - I3udu -{(3- ^)vd^, 

X = X2 = e^\a^du + a^d^ + (3udu + (33) 

X = X3 = e^'+P-^i^a^du + a^^^) 

where the Greek letters again denote arbitrary parameters and p ■ x = p^x^. 

The next step is to specify all non-equivalent sets of arbitrary constants in 
using the equivalence transformations fll2l) . 

If the coefficient for udu (or vd^) is non-zero then translating u (or v) we 
reduce to zero the related coefficient (a^) in X^ and X^; then scaling u (f ) 
we can reduce to ±1 all non-zero a" in (1331) . In addition, all operators (l33ll are 
defined up to constant multipliers. Using these simple arguments we come to the 
following non-equivalent versions of operators 



(34) 



Xf ^ = 2pD - udu + (/i - l)f5^, 
Xf ^ =2D + vd^ + vdu, Xf ) =2D- ud^ - d^, 

To describe two-dimensional algebras A we represent one of the related basis 
element X in the general form flH^ and calculate the commutators 



Y = [p°, X] - 2pP°, z = [p°, r], = [X, r] 

where is operator given in ( |26l) . After simple calculations we obtain 

r = Fi(M9„ + t;9,) + 51 9^ + 525^, Z = Fi4(Ma„ + t;5,) + 5ia„ + 52a,, 

= 2ptZ + px,{Bl^du + Bl^d,). ^^^^ 
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By definition, Y , Z and W belong to A. Let 7^ then we obtain from fl35p : 



/x^O: B1, = Fu = Bl = Q, (36) 
/i = 0: Fu = aFt + rB^, B^, = ^'^ Ft + f3^' B',. (37) 

Starting witli (IHUj) we conclude that up to translations of t the coefficients F 
and -B" have the following form 

F = at or F = (3; 5" = + a" if /i ^ 0. 

If F = crt then the change 

Ua ^ M.e-'^* - —t (38) 

reduces the related operator ( l22l) to of ( l33i) for = 0. 

The choice F = j3 corresponds to the following operator ( l32l) 



X = = - 2t(a^(9„ + a^S^) (39) 



where X^ is given in (!33|) . 

Thus if one of basis elements of two dimension algebra A is of general form 
with /i 7^ then it can be reduced to X^ with /? = or to generator (IH^ . 
We denote such basis element as e^. Without loss of generality the second basis 
element of ^ is a linear combination of operators Xg'^'' and Xg'^'* (1341) . Going 
over possible pairs (e^,e^) and requiring [e^,e^] = a^e^ + a^e^ we come to the 
following two dimensional algebras: 



Ai =< 2D + vd^, X^°) >, A2 =< Xf \ Xf ^ >, 

A3 =< X[^\X^^^ >, A4 =< X[^\ Xg^-* >, 

=< X(\), Xf ^ >, Ae=<2D + 2vdy + udu + vtd^, xf^ >, 
=< 2D + 2udu + 3t;9^ + ?,utd^, X|\) > . 



(40) 



The form of basis elements in (1401) is defined up to transformations (]T2l) . (l38i) . 

If A does not include operators (1321) with non-trivial parameters 11 then in 
accordance with flHSl) its elements are of the following form 



= (na„ + t;9,) + Bl^d^ + a = 1, 2 (41) 

where F*^"^ and -B^^^^ are solutions of (1371) . 

Formulae fjlUl) . fBTl) define all non-equivalent two-dimensional algebras ^ 
which have to be considered as possible symmetries of equations (El). We will 
see that asking for invariance of ([H]) w.r.t. these algebras the related arbitrary 
functions are defined up to arbitrary constants, and it is impossible to make 
further specification of these functions by extending algebra A. 
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6 Group classification of equations ([81) 



We suppose that parameter p in (IHD be nonzero. Then, scahng independent 
variables t, x we can reduce it to p = 1. 

To classify equations ([8]) which admit one- and two- dimension extensions of 
the basis invariance algebra (1261) it is sufficient to solve determining equations 
f l24p for with known coefficient functions -B" and F of symmetries fl32p . These 
functions are easily found comparing ( l22l) with (HOl) and (HTl) . 

Let us present an example of such calculation which corresponds to algebra 
Ai whose basis elements are = 2tdt + Xadx^ + vd^ and x'f''' = udu + f <9^,, refer 
to ( HOl) . Operators and Xg""* generates the following equation (l24l) . ( l25l) : 

/' = -«/.^ f = -^«/.' (42) 

and 

r = (n9„ + t;9,)r, a = 1,2. (43) 

General solution of (03]) is: = uF^ (^) , p = uF^ (^) were F^ and 
are arbitrary functions of ^. Solving ( l42l) for such functions and we obtain 

/I = au\-\ f = \u\-\ (44) 

Thus equation ([H]) admits symmetries Xq^"* and X^ provided and are 
functions given in (jS]). These symmetries are defined up to arbitrary constants 
a and A. If one of these constants is nonzero, than it can be reduced to +1 or 
— 1 by scaling independent variables. 

In analogous way we solve equations ^IM corresponding to other symmetries 
presented in (1341) and ( l40l) . At that we do not consider functions and which 
are ether linear in m, v or correspond to decoupled systems (IHl) (i.e., when 
and p depend only on u and v correspondingly). The classification results are 
presented in Table 1. 

In the fourth column of the table symmetries of the related equation ([H]), 
are presented together with the additional equivalence transformations (AET) 
which are listed in formula (|T5|) : the numbers of AET from the list (1151) are 
given in square brackets. Greek letters denote arbitrary real parameters which in 
particular can be equal to zero. Moreover, without loss of generality we restrict 
ourselves to ?7 = 0, 1, 5 = 0, ±1, e = ±1. 

In Table 1 D is the dilatation operator given in (ITSl) . x = {xi,X2, ■ ■ ■ ,Xm-i), 
\E'(a;) is an arbitrary function of spatial variables; \E'^(x, Xm — f]t) ^m(^' ^) 
solutions of the Laplace and linear heat equations: 

A„^^ = /i^^, (I - A„_i)<l>^ = 
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Table 1. Non-linearities and symmetries for equations (8), (9) with 

p=l 



No 


Non-linearities 


Arguments 
of 


Symmetries 

and AET Eq. (15) 


1 


fl = ^2.+l^l ^ 
r'} 1J-\-^ 7-i9 




2vD — udu + — '^)vdy 
[AET 1 if z/ = 0] 


2 


f = F^v-\ 

£1 r-<2 —1 
= F^f ^ 


M — In f 


2D + vdy + ridu 


3 


/I = + elnM), 
= + e Inw) 


V 

u 


e*^* {udu + "i^c^ii) 


4 


/I = u'F\ 


V — Inu 


2D — udu — du 


5 


p = F' + ^,u, 

x2 7^2 1 

f^ = F^ + TjU 


V 


e~'i^"-^^{t,x)du 

[AET 2, 22 if 7] = n = 0\ 


6 


f'^ = F'^ + vv 


u + rjv 


e''^-^i,{x,Xm-vt){dv -vdu), 
[AET 11 if u = jj, = 0\ 


7 


rl 2—1 
/ = jlU V ^ 




2D + vdy, udu + vdu 
[AET 1 ] 


8 


fl ^ ^g-2n^ 
x2 —11 

f^ = r]e " 




2D + vdu + du, '^{x)du 
[AET 3] 


9 


x2 r 2i) 




2D - udu - du, $o(t, x)du 
[AET 2, 22] 


10 






2uD - udu + l)vd^, '^{x)du 
[AET 3] 


11 


xl Sv— 2 




2(z/ — 1)D — uudu — vd^, 'Po[t, x)du 
[AET 2, 22; ] 


12 


/i = ^, /2 = lnt. 




2D + 2vd^ + udu + td^, '^{x)d^ 
[AET 3; & 7 it u = 0] 


13 


/-"^ = Inf, = i/f 5 




2D + 2udu + 3vdy + 3tdu, 
%{t,x)du 

[AET 2, 22; & 8 if u = 0] 



7 Group classification of equations ([2]) 

In this Section we present the classification results for coupled systems of equa- 
tions ([2]). The related classifying equations are given by relations (!20l) . 

Like in Section 5 we first describe all non-equivalent low dimension algebras 
of the main symmetries for equation ([2]). Non-equivalent realizations of these 
algebras (together with detailed calculations) are present in the Appendix. Using 
found realizations of algebras A and solving the related classifying equations 
we easily complete the group classification of equations ([2]). 
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We will not reproduce here the related routine calculations but present the 
results of group classification in Tables 2-10. Besides symmetries and the related 
non-linearities, the additional equivalence transformations which are admissible 
by particular classes of equations ([2]) are indicated there. The symbols G^, G'^ 
and K denote generators listed in (ITSl) . V'm = ^^(x) are arbitrary 

solutions of the linear heat equations and Laplace equation: 

and \E'(x) is an arbitrary function of x. The Greek letters denote arbitrary pa- 
rameters. Moreover, up to equivalence transformations we restrict ourselves to 
e = ±1, = 0, 1 and 5 = 0, ±1 . 



Table 2. Non-linearities with arbitrary functions and symmetries for 
equations (2) with arbitrary a ^ 



No 


Nonlinear terms 


Arguments 
of F\ F^ 


Symmetries 


AET 

Eq.(15) 


1 


p = u''-^'F^ 




uD — udu + fJ'Vdy 
for any /i, u, 
& Ga for z/ = 0, 
a/i = 1 


14, 

p = jlOJ 

±f V = y) 


2 


p = u{F^ + e\nu), 
p = viF^ + e\nv) 




e^* (udu - fxvd^) 
for any /i, & 

Ga if a/i = —1 




3 


P = v''F\ 
p = v^+^F^ 


u — \nv 


uD — vdy — du 


4,/i = -l 
if z/ = 


4 


p = F^ + eu, 
P = F'^v + euv 


u — \nv 


e^* (vd^ + du) 




5 


P = 0, P = F\ 
a ^ 1 


u 


D + vdu, ipodv 


3 


6 


= F\ 
p = F^ + 6v 


u 


ijjsdu 


3 if 5 = 0; 
6 if a = 1, 

F^ = 6u 


7 


p = F^ + 6u, 
p = F^ + av, 
a ^ 1 


V — u 


e^'^H^ix) {du + d,) , 

(l-a)' 

K = a + afi 




8 


P = e"F\ 

P = e"F2 

1] = if a = 1 


V — rju, 


D-du- r]dv 





In Table 3 A denotes the characteristic determinant, A = \{fJ^ — i^Y + Xcr. 
Additional equivalence transformations are specified in the third column. 
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Table 3. Symmetries of equations (2) with cirbitrary a 7^ and 
non-lineeirities p = u {iilnu + Xlnv) , p = v {ulnv + alnu) 



No 


Conditions 


Symmetries and AET Eq.(15) 


Additional symmetries 


1 


A = 0, (7 = £, 

Li — ly 


e^^vdv, e'** {udu + stvdu) 
[aft 90 1 f // — nl 


none 


2 


\ = 0,a = e 
/i ^ z/, (a - If 
W 


e^* ((/i — z/) -uc^u + (jfflu) , 
e'^*t;9^ [AET 14,tu = 
if /XI/ = 0] 


Ga if i/ = —aa, /x = 0; 
ii II ^ 0, II — ly = aa 


3 


A = 0,A(7^0, 
A2 + ^2 = 1 

^ + z/ = 20 


= e''\2\udu + {v - IJi)vdu), 
2e^^vdu + tX2 [AET 14, 
acu = —up if /i + z/ = 0] 


Ga if // = — i/, \ = av 

Ga'^iv^ 

2A = a(z/ — /i) 


A 


A = l, 
u;± = ± 1 


[AET 14, era; = -z/p 
if /xz/ = Act] 


It it 7 / / / — A /T A — — n I / 

VJTq; 11 /\\J ^ /\ XAjUj 

Ga if py 7^ Act, 
A = a{y — 11 + aa); 


5 


A = -1 


e^^^ {2X cos tudu+ 

{{v — jj) cos i — 2 sin t) vdy) , 

e^* (2Asin tudu+ 

( (z/ — /i) sin t + 2 cos t) vS^) 


none 



Table 4. Non-linearities with arbitrary parameters and extendible 
symmetries for equation (2) with any a 



No 


Nonlinear terms 


Main 

symmetries 


Additional 
symmetries 


AET 
Eq. (15) 


1 




I'D — tic}„ 


Ga if au = —n, 
& Xif 
i/m(l — a) = 4; 


14, 

ucu + lip 
= 


V'o^u if CT = 0, 
u = -l, & 
if /X = a, & 
if a = 1 + f ; 


2; 14, 
uj = pp 


2 


z/2 + (a - 1)2 ^ 


uD — udu ~ l^vdy. 


Ga if z/ = 0, 
a/x = 1 


3; & 14, 
p = poj 
if z/ = 


3 


P = lnu, a ^ 1 


D + udu + vdu+tdy, 

ijjodv 


+ tc?^ 
if (5 = 


3,7,9 


4 


P = Sulnu, 
p — j/v + In 14, 

a ^ 1, 

z/2 + 5 V 


ipudv 


e-"' {udu + tdy) 
iiv = 5 




e'^ {{5 - v)udu 
+dy) a 


if (5 = 


5 


y2 _ 


vD - vdu - du. 


{u — 5t)dy if 
u = 0,a = 1 


3, & 9, 4, 

yU. = — 1 

if i/ = 
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Classification results present in Tables 4 and 5 are related to systems with 
arbitrary values of a presented by equation (1141) (if not specified in the second 
columns of the tables). 

The non-linearities given in Table 4 are defined up to arbitrary parameters. 
For some values of these parameters the related equation ([2]) admits extended 
symmetries indicated in Column 4 of the table. 

In Table 5 non-linearities for equation are classified whose symmetries are 
fixed for all admitted values of parameters. 



Table 5. Non-linearities with arbitrary parameters and 
non-extendible symmetries for equations (2) with arbitrary a 



No 


Nonlinear terms 


Symmetries 


AET Eq. (15) 


1 


J — [u-\- V) , 

f^ = fi{u + vY^^ , a ^ 1 


Till HI /I ! 

iJU — UjUu — u(J^, 


11,77=1 


2 


/ - e , 

/2 = ee^ a ^ 


D-du, ipodu 


2 


3 


/I = 5e"+^ 

/ - ae ^ , a 7^ 1 


D-du, 


11, ^ = 1 


4 


a ^ 0, + /i V 


D - du, vd^ - fidu 


4, if 
a = 


5 


P = ee\ p = u 


D + vdy — du — tdy. 


3 


6 


P = e\n{u + v) , 
P = ulnlu + v) , 


^o(a:) {du - d,) , 
e{a -1){D + udu + vd^) 
+ {{a + ev)t 
+ {du - du) 


11, 77=1 


7 


p = eu^+\ /2 = lnu, 


V {D + vdu) — udu — tdu, 
Adv 


3, & 7 
if = 


8 


P = (/i — In M + uv, 
P = — z/^ Inw + (yU + u)v 


X3 = e'^* {udu + i^du) , 
tXs + e^'du 


5, K = —V 
if /i = 


9 


P = {fi — i')u\iau + uv, 
/2 = (1 _ u^)lnu 
+ {fx + u)v 


Xf = e^f'^^'>\udu 
+{^±l)du) 


5, K = /i — z/ 

if /i = ±1 


10 


P = {jj, — u)u In M + uv, 

/2 = (/i + V)V 

- (1 + Inw 


e^* {cost{udu + I'dy) 
— sin tdu) , 
e^* {smt{udu + vdu) 
+ cos tdu) 





In Tables 6-9 and 10 the additional symmetries are presented which corre- 
spond to the specific values a = 1 and a = of the diffusion coefficient. We use 
the following notations here: R = yu? + v^, z = tan~^ ^. 

In Table 6 the AET are given in square brackets and placed in the last column. 
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Table 6. Additional non-linecirities with cirbitrciry functions and 
symmetries for equations (2) with a — 1 



No 


Nonlinear terms 


Arguments 
of F^, F2, F 


Symmetries 

and AET Eq. (15) 


1 


= uF^ + 6riv, 
+uF^ + vF^, 

u // y- u 




e"^* {udy + ri{udu -\-vd^)) 

for any r^^Fi 
& "ipadv for ?7 = 0, 

= a - 5 ^ 0, [AET 3 
if F^ = -(5, 77 = 0] 


9 


P = u^+^F^ 

P = u"" {F^v + F^u) 


-n- 


rji^vD — udu — vdy) — udy 
[AET 6 if 77 = 0] 


3 


p = uF^ + vF'^ 

+ez [jiu — v) , 
p = vF^ - uF'^ 
+ez {ixv + u) 




e^* (liRdR + d,) 


4 


p _ QVupi^ 

p = e^"(F2 + F^u) 


2v-v? 


r]D - udy - du 
[AET 13 if 77 = 0] 


5 


P = eu + F\ 

P = eu^ + F^u + F2 






6 


P = Fu, 
P^Fv 


u 


[AET 6; & 14,a; = 0] 


7 


P ^rj, p^Sv + F, 


u 


iisdy, e^^{u-r]t)dy 

for any 7/, 5, & 
D + 7;9.y for 1] = 6 = 
[AET 3 if 5 = 
& 6 if 77 = 0] 


8 


fl ^ gAz / pi I p2 \ 

/2 = {F^v - F^m) 




\D + V {udu + vdy) 
-udy + vdu 

[AET 15, (7 = 1 if A = 0] 


9. 


/I = elF^u, 

p = el {Fh, + F2) 


u 


D — ud-v 


10 


/I = p = (u + ()> + F 


u 


e^^udy, e^^{dy + tudy) 


11 


p = (u + iy)v + F 


u 


e('^-i)* {ud, - d,) 


12 


p = {u' + l), 

p = {u + v)v + F 


u 


e^* {costudy — sintdy) , 
e"^ {smtudy + costc?^) 



In the following Tables 7 and 8 A = \{jjL—vY-\-\a. Symmetries and additional 
equivalence transformations are specified in the third column; AET are given in 
square brackets. In the last columns additional symmetries are specified which 
are valid for some particular values of parameters defining non-linearities. 
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Table 7. Symmetries of equations (2) with a — 1 and non-lineEirities 

— {/lu — av) In + z{Xu — uv) , p — {/iv + au) lnR + z{Xv + vu) 



No 


Conditions 
for coefficients 


Symmetries and AET Eq. (15) 


Additional 
symmetries 


1 


A = 0, 

jj, ~ V — 5 


e''{RdR + atd,) 
[AET 16 if n = 0] 


G« if (7 = 0, ^ 7^ 


2 


A = 0, // ^ u, 
li^ + u^ = l 


e-^'d,, e^' {ad, + (ii - u) ROr) 
[AET 15 if fiu = 0] 


Gq if = (T = 0, 
G« if /X 7^ 0, (J = 


3 


A = 0, 

X = e, 

n + u = 2n 


X5 = e^" {26RdR+{i^-li)d,), 

2e^'d, + tXs 

[AET 15 if i/ + = 0, 

& 1, 17 if /X = 1/ = 0] 


Gq if /X = 1/ = 0, 
Gq; ii /I — f ^ 


4 


X^O, A = 1 

u}± = n±i 


e-+*(Ai?a^+(a;+-/x)9,), 
e--*(Ai?9H + (cu_-/i)a,) 
[AET 15 if = Act, 
& 1 if // = (7 = 0] 


Ga ii cr — II — 0, 
Ga if CT = 0, /i 7^ 


5 


A = -1 


exp(nt) [2 X cos tRdR 

+ ((z/ — /i) cos t — 2 sin dz] , 

exp(fit) [2Asinti?9R 

+ ((z/ — /i) sin t + 2 cos t) dz] 


none 



Table 8. Symmetries of equations (2) with a = 1 and non-linearities 

f = Xv + iiulnu, f = X^ + {au + ^iv) hiu + uv, A^ + 7^ 



No 


Conditions for 
coefficients 


Symmetries and AET Eq. (15) 


Additional 

symmetries 


1 


A = 0, 
11 = V = 5 


e^^ud^, e^\RdR + atudy) 
[AET 10 if /i = 0] 


if 11 = 


2 


A = 0, // ^ V, 

//2 + = 1 


e^* ((/i — u)RdB + audy) , e'^^udv 
[AET 18 if //// = 0] 


V'„dr if /( = 0. 


3 


A = 0, 

A = 

/X + z/ = 2f] 


= e^^* (2ei?(9R + {iy- ii)ud^) , 
2e^*?7,a^ + tX^ [AET 18 if 
ji = -V, k 1, 19 if /i = z/ = 0] 


D + udu, Ga 
if = z/ = 0; 
G'^ if ii = iy^O 


4 


A 7^0, 
A-1, 

a;± = ± 1 


e"+* (Aiia.^ + (a;+ - ii)ud,) , 
e'"-* (Ai^^ij + (a;_ - //)ii9^) 
[AET 18 if 111/ = Act, 
& 1 if /i = a = 0] 


Gaii a = n = 0, 

Ga if (7 = 0, /X 7^ 


5 


A = -l, 


e^"[2Acosti?9R 

+ ((z/ — /x) cost — 2 sint)M(?^], 

e^*[2Asinti?9fl 

+((1/ — ijl) sint + 2cost)it9„] 


none 
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Symmetries presented in Table 9 are related to equations ([T]) with the unit 
diffusion matrix only. 



Table 9. Additional non-linearities with arbitrary parameters and 
symmetries for equation (2) with a = 1 



No 


Nonlinear terms 


Symmetries 


AET Eq. (15) 


1. 


/I = s, p = «^ 

6 = if u = 2 


D + udu + (z^ + l)vdjj, 
-ijjody, [u - 6t)dy, 
for any z/, 6, & 
udu + i^vdv for 5 = 
& 9u + 2tudy for z/ = 2 


3,9;& 14, 
p = uuj if 
6 = 0, k 
21 if z/ = 2 


2. 


/I = eu, p = u\ 
z/ ^ 0, 1 


udu + i^vdv, ipody, e~'^^udu 

for any z/, & 

d^^{udu + edu) f or z/ = 2 


3;& 

14, p = uui 






D — vdv, udu + vdy, Ga, 




3. 


P = ^v, P = -'i 


K + {2 - 7]){t{rjudu 
-(2 + r])vdu) - udy) 


1 


4. 


p = eu^+\ 

p = eu^v, z/ 7^ 


uD — udu, vdy, udv 
for any z/, & (1 + tu)dv 
for u = 1 


Q;k 
14,u; = 




p = Su''+\ 






5. 


nO 1/ / r , /-r\ 

= u" {5v + pu") , 
+ (T 7^ 0, juz/ 7^ 


uD — udu — crvdu, ud^ 


6 


6. 


P = 5 {2v - u'Y^'^ , 

P = 6u {2v - u^y-^-' 
_L,, (2v - 

\ LI LI J 


2z/iJ' — udu — 2vdy, 
du + udy for any z/, fi 
k 2t{udu + du) + du 
for z/ = 2, /i = 


13; & 14, 
p = 2uj 
if z/ = 


7. 


£1 r n 

/ = cie", 
J — ue 


6{D - du)-udu, ipody 


3 


8. 


p = ne^^-^\ 

p = (nu + p) e^''-"' 


2D — dy, du + udu 


13 


9. 




D — udu, uD — udu — vdu 


12 


P = e^[Sv + au)u'' 


for any v, & Ga for z/ = 


10. 


p = e'''R''{Su- fxv), 
p = e'''R''{5v + fiu) 


aD — udu -~ vdu, 

uD — udu + vdu for any a, 

kGa for a = 


15 


11. 


P = eu \nu, 


e*^* (m^u + vdu) , vdu, vdu. 


6; & 


P = evlnu 




14,cu = 


12. 


P = 5, 
P = \nu 


D + m9„ + + tdu, i>odu 
{u — 6t)du for any 6, 
& udu + ^f^j; for 5 = 


3,9, 
& 6, 7 
if 5 = 
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Table 9. Continued 



13. 


P = eu^{v — Inu), 
/i 7^ U 


jlD — UOu — Oy, UOy 

for any /i 7^ 0, 

&c + LUOy lor /i — i 


6 


1 A 
14. 


/I = 51n(2w-u2), 

J — (J\ZV — U ) ' 

+6u ln(2f — -u^) 


+ udu + 2f 9^ 

+101 [Ou+UOv ) , 

du+udy 


id 


15. 


f2 _ u+l ] 

/ (X 111 (X 


uD — {udu + vdu + £udy), 

(foOy 


3 


16. 


/I = £u''+\ V ^ 1 

P = eu^v + ulnu 


uD — udu — tudu 
—(1 — iy)vdy, udy 


6; & 5, K = e 
if z/ = 


17. 


P = 2v - u^, 

p = {fi + u)\2v-u^) 

2 

— ^u, /X 7^ 


= e'^* {2du + 2udu 
+fxdu), tX' + e^^'du 




io. 


P = 2v- u^, 

J — [fJ' + U) [ZV U ) 


X± = e(M±i)t (25„ 

+2udu + {l2±l)du) 


13 if 


19. 


P = 2v- u^, 

p= 


e^^ {2cost{du + udu) 
+ cost — smt)du) , 
e^^* {2smt{du + udu) 
+ sint + cost)du) 





We did not consider decoupled systems ([2]) whose symmetries can be easily 
found using the classification results of Dorodnitsyn p!8] for a single diffusion 
equation. We also did not specify the case of hnear systems ([2]) when 

P = vu + fiv + a, p = an + Xv + uj. (45) 



Equivalence transformations ffT^ and 1-3 of ffT^ make it possible to specify 
values of parameters in fH3|) by imposing the following conditions: 

a; = ti; = A = 0; fia = or fi = icr. (46) 

Moreover, if the diffusion matrix A is proportional to the unit matrix then equa- 
tion ([2]), fHHl) can be reduced to the case P = P = 0. 

The classification results present in the tables are valid also for equations 
([2]) whose r.h.s. have the form ( H5i) . ( H6l) . However, to save a room we did not 
indicated the standard additional symmetries of linear equations, i.e., Udu and 
Vdu where U and V satisfy the relations 

Ut-/^U = p, Vt-aAV = p. 

The following last table completes the classification results for the case of 
singular diffusion matrix. 
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Table 10 

Additional non-linearities and symmetries for equations (2) with a = 



No 


Nonlinear terms 


Argu- 
ments 
of F\ F^ 


Symmetries and AET (15) 
[in square brackets] 


1 




V — u 


[AET 11, ?7 = 1 
if ;x = 5 = 0] 


2 


F2 = 1 if r/ = 


V — rju 


D - du - r]9v 


3 


p1 7-t1 

= F\ 
f^^F^ + r]v 


u 


[AET 3 if ?7 = 0] 


4 


P = vF\ p = F' 


u 


D + vd^ 


5 


p = uF^ + 5uv, 
p = F^ + 5v 


V — \\iu 


[AET 5 if 5 = 0] 


6 


/> "1 t J_l_ 1 T — I 1 

p = u^-^^F^, 

p = u^F^, u^O 


V — \\iu 


uD — udu — 


7 


P = F^ + uu, p = ri 


V 




8 


p = u^+^F\ p = 


V 


vD — udu f 01" any v 
& iIjqOu if = — 1 
[AET 2, & 14, p = 
if z/ = 0], 


9 


p = v^+\ p = 6 




D + vdu- Xudu, ipodu 
[AET 2] 


10 


p = 5e", p = 




D - du, ^{x)du 
[3; & 4 if (5 = 0] 


11 


p = \nv, p = e 




D + udu + vdy + tdu, 
^odu [AET 2] 


12 


p = Su''+^v~\ 
P = u'' 




D + vdv, I'D — udu 
[AET 14, uio = p 



In Table 10 is an arbitrary function of v. In addition to the equivalence 
transformations indicated in the fourth column, all the corresponding equations 
([2]) admit the AET m — > m, w — > (f{v) where if is an arbitrary function of v. 



8 Discussion 

We have carried out the group classification of systems of coupled reaction- 
diffusion equations ([2]) with a diagonal diffusion matrix. The classification results 
are present in Tables 2-10. Moreover, symmetries of equation ([8]) with a singular 
diffusion matrix and additional first derivative terms are presented in Table 1. 
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The list of non-equivalent systems ([2]) appears to be rather extended, espe- 
cially for the unit diffusion matrix. Equations with invertible and non-unit 
diffusion matrix A have an essentially shorter list of different symmetries. If 
the diffusion matrix is singular the number of inequivalent equations appears 
to be the smallest one which is caused by the powerful equivalence relations 
u u,v —>■ (f){v) where (/> is an arbitrary function of v. 

More exactly, if matrix matrix A be of type 1, equation ffT^ . then there exist 
9 non-equivalent classes of equations ^ defined up to arbitrary functions and 19 
classes of such equations defined up to parameters. The related non-linearities 
and symmetries are presented in Tables 4,5 and 10. The presented extensions of 
the basic symmetries (126|) have dimensions from 1 up to 3 and include neither 
Galilei generators Ga nor conformal generators K. 

In addition, in Table 1 thirteen classes of equations with a singular diffusion 
matrix and first derivative terms are presented. 

For the case when matrix A is of type 2, equation ([HI), we indicate in Tables 
2-5 ten classes of equations defined up to arbitrary functions and thirty five classes 
of equations defined up to arbitrary or fixed parameters. Among them there are 7 
Galilei invariant systems, whose r.h.s terms are given in Table 2, Item 2; Table 3, 
Items 2,3,4 and Table 4, Items 1,2. In addition, there exist two systems of type 
([2]) with a diagonal (but not unit) diffusion matrix, which are invariant w.r.t. 
extended Galilei algebra spanned on P^, J^^ ( l26l) dilatation operator and also 
generators Gq, K (jTSIl . These equations correspond to the non-linearities present 
in Table 4, Item 1 and have the following form 

4 

ut- Au = Xu (uf"") , 

4 

Vt - aAv = av [uv^'') ""'^-"^ 

and 

ut — Au = Xv^^ , 

Vt - (l + — I Av = 0. 
\ m J 

Finally, if the diffusion matrix is the unit one then we indicate 98 non- 
equivalent classes of equations, among them 21 including arbitrary functions and 
14 admitting Galilei generators. There is the only equation admitting extended 
Galilei algebra, the related non-linearities are given in Table 9, Item 3. 

Consider examples of well known reaction diffusion equations which appear 
to be particular subjects of our analysis. 

The CGL equation with /3 = can be rewritten as 

Ut - A2U = u+ {u"^ + v'^){av - u), , . 

Vt — A2V = V — {u"^ + v'^){v + au). 

where u and v are real and imaginary components of the complex function W. 



21 



The r.h.s. of equations (1471) has the form presented in Item 8 of Table 6 (with 
A = z/ = 0), and so in addition to basic symmetries < do, di, d2,Xid2 — X2di > 
this system admits the symmetry 

X = ud^- vdu. (48) 

Using the anzatse 

^ _ gj(a;i cose+X2 sine)^^ ^ _ gi(a;i cos e+X2 sinfi) ~ 

where u and v are functions of t and uj, uj = Xi sin 9 — X2 cos ^, ^ is a parameter, 
the system f H7|) can be reduced to the form 

u - u^ui = (u^ + v'^)iav - ii), vt~ Vuuj = {u^ + v^){v + au). (49) 

Main symmetries of the reduced equation ( H9l) appear to be more extended 
then of the CGL one. As is indicated in Item 11 of Table 9 equation (H9!) admits 
symmetry (HSj) and also the following one: 

X2 = 2D — udu — vdy. 

The primitive predator-prey system (jll) is a particular case of equation ([T]) 
with the non-linearities given in the first line of Table 2 where however — /i = u = 
1, = —F"^ = ^ . In addition to the basic symmetries < dt, dx > this equation 
admits the (main) symmetry: 

X = D — ud„ — vdy. 



The X — u react ion- diffusion system ([5]) and its symmetries was studied in 
paper [15]. Our investigations confirm and complete the results of [15]. First we 
recognize that this system is a particular case of ([T]) with non-linearities given 
in Item 11 of Table 6 with fi = u = 0. Hence it admits the five dimensional 
Lie algebra generated by basic symmetries (1261) with /i, i/ = 1,2 and also the 
symmetry (1481) . This is in accordance with results of paper [12] for arbitrary 
functions A and lo. Moreover, using Table 9, Item 11 we find that for the cases 
when 

X{R) = XR", 00 = aR" (50) 

equation ^ admits additional symmetry with respect to scaling transformations 
generated by the operator: 

X = vD — udu — vdy. (51) 

The other extensions of the basic symmetries correspond to the case when X{R) = 
fi\n{R),u{R) = a\ia{R), the related additional symmetries are given in Table 7, 
Items 1, 2, 5 where u = X = 0. 
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Consider now the system (jH]). This system admits the equivalence transfor- 
mation 1 fllSI) for p = —uo. Choosing p = 2k we transform equation to the 
form ([T]) where a = —1, = —2u'^v and = 2v'^u. The symmetries corre- 
sponding to these non-hnearities are given in the first hne of Table 4. For m = 2 
the symmetries are the most extended and include two dilatations, two Galilei 
generators Ga a = 1, 2 and the conformal generator K. All these symmetries 
except K are valid for other numbers m of independent variables. 

Symmetries of equations ^ for m = 1 were investigated in paper [TT] whose 
results are in accordance with our analysis. 

The results of the present paper related to non-degenerated diffusion matrix 
can be compared with those of [1| and ^Sj, ^Qj. 

Paper ^ was apparently the first work were the problem of group classification 
of equations with a diagonal diffusion matrix was formulated and partially 
solved. However the classification results presented in include only a small 
part of ones presented in Tables 2-10. 

In papers [8], [9] Lie symmetries of the same equations and also of systems 
of diffusion equations with the unit diffusion matrix were classified. The results 
present in those papers are much more advanced then the pioneer Danilov ones, 
nevertheless they are also incomplete. In particular, the cases presented above 
in Items 13-15 of Table 6 and Items 1, 2 of Table 7 were not indicated in [9j. 
Moreover, many of equations treated in [9] as non-equivalent ones, in fact are 
equivalent. For instance, all versions 14, 15, 18 and 20 from Table 4 present in 
[H] are equivalent one to another. 

Notice that the results related to the group classification of systems of non- 
linear systems of reaction-diffusion equations are presented in very compressed 
form and discussed in the survey [19] . The principally new points of the present 
paper in comparison with [12] are the following ones: 

• In the present paper we give the completed list of admissible equivalence 
transformations f|T5l) for all classified equations ([2]) whereas in [19j only an 
a priori fixed subclass of equivalence transformations was discussed. 

• For any particular system of equations ([2]) whose nonlinear terms are given 
in the classification tables the admissible equivalence transformations are 
specified and presented explicitly at the same tables while in [T^ the general 
(incomplete) list of such transformations was presented only. 

• We use our knowledge of all admissible equivalence transformation to reduce 
the number of non-equivalent versions of systems ([2]) to absolute minimum. 
In particular many of quantities which define nonlinearities and are treated 
in [19] as arbitrary parameters are reduced to S = 0, ±1, e = ±1 or = 0, 1 
and possible values of parameter a in the diffusion matrix A are reduced to 
ones given by equation (fT^ . . 
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• Summarizing, in the present paper the problem of group classification of 
systems of reaction-diffusion equations is solved completely whereas all 
previous publications jl]-[9] and [19] can be treated only as steps to the 
complete solution. 

Thus we present group classification of reaction-diffusion systems with a di- 
agonal diffusion matrix. Such systems with the square and triangular diffusion 
matrix has been classified in paper [Ij and preprint [3J respectively. The results 
of papers [1], [3] and the present one consist in the completed group classification 
of systems of two coupled diffusion equations with the general diffusion matrix. 



A Appendix. Algebras of main symmetries 

Following pp we first specify all non-equivalent terms 

N = C'^'ubd^^ + B-d^^. (52) 

where summation from 1 to 2 is imposed over the repeated indices and we again 
use the notations Ui = u,U2 = v. 

Let fl3^ be a basis element of a one-dimensional invariance algebra A then 
commutators of N with P° and should be equal to a linear combination of 
and operators ( l26l) . This condition presents three the following possibilities [I]: 

1. = B" = fi", 

2. C"^ = e^*/i"^ = e^V", (53) 

3. C"^ = 0, P" = e^*+^-^>" 

where A, and uj are constants. 

Like in [T] to classify all non-equivalent symmetries 0531) we use their isomor- 
phism with 3x3 matrices of the following form 

\ 

fi' fx'^ . (54) 

Equations ([2]) admit equivalence transformations ( |T2l) . The corresponding 
transformation for matrix (ISH) are 



1 

g' = UgU-\ U = \ K^^ K^^ I (55) 
were K"-^ are the same parameters as in (1121) . fll3p . 
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For the case of equation with a 7^ 1 matrices fi and K in (13^ . (1331) are 
diagonal, and up to equivalence there exist there exist three matrices (13^ . namely 






g' = \ 1 \, (7=100, g' = \ \ \. (56) 



In accordance with flTTj) . fl3^ . fl3Hl) the related symmetry operator can be rep- 
resented in one of the following forms 

^(fc) = f^D ~ 2{g'')hcUcdu^, X^,) = e^\g'')bcUcdu„ 

(57) 

^3 ^ gAt+...x (5^^ + ^5„J , A; = 1, 2, 3. 

Here (5^^)60 are elements of matrices (!3U|) . 6, c = 0, 1, 2, m = column (1,mi,'U2)- 

Formulae ( |571l and ( !56l) give the principal description of one-dimension alge- 
bras A for equation ([2]) with a 7^ 1. 

To describe two-, three- and four- dimension algebras A we first classify the 
corresponding algebras An^s of matrices g fl3^ where index n indicates the dimen- 
sion of the algebra and s is used mark different algebras of the same dimension n. 
Choosing a basis element of A2^s in one of the forms given in fl3U|) we find that up 
to equivalence transformations (l54l) there exist six two dimension algebras with 
basis elements < 61,62 >: 



^2,1 : ei = fifjg), 62 = g^; ^2,2 : 61 = fiffp), 62 = gf^y, 
^2,3 : ei = ^^ 62 = ^^0), 



^2,4 : ei = g\ 62 = 5-^ ^2,5 : 61 = 62 = g 
^2,6 : ei = 63 = 5(1*0) 



(58) 
(59) 



where g^^^ = g^\x=o, g^^-^ = g^\\=i, gf^) = 9^\x=o, and 

/000\ /000\ 
= , / = 1 . (60) 
\001/ \000/ 

Algebras ( !58l) are Abelian while algebras (|59l) are characterized by the follow- 
ing commutation relations: 

[61,62] = 62. (61) 



Using fl58l) . fl59|) and applying arguments analogous to those which follow 
equations flH^ we find pairs of operators flTTj) forming Lie algebras. Denoting 

d 

Cq = (6Q)a6Ub7T::-, = 1,2 
dUa 
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we represent them as follows: 

< HD + ei + 1^^2,62 >, < fiD + 62 + 1^^1,61 >, , s 

<HD- ei, uD - 62 >, < F^ei + G^es, F^ei + G^e2 > ^ ^ 

for ci, 62 belonging to algebras and 

< /xD - ei, 62 >, < fiD + ei + i^te2, 62 > (63) 



for ci, 62 belonging to algebra 

Here {F^, G^} and {F^, G^} are fundamental solutions of the following system 

Ft = XF + vG, Gt = cTF + 7G (64) 

with arbitrary parameters A, i/, o", 7. 

The list (1621)- flU^ does not includes two-dimension algebras whose basis is 
< Fe„, Ge„ > (with F, G satisfying ([MD) or < /iL) + Ae''*+'^-^e<„ e^*+'^-^e„ > which 
are incompatible with classifying equations ( l20l) . In the following we ignore all 
algebras A which include such subalgebras. 

Up to equivalence there exist three realizations of three-dimension algebras of 
matrices CT. (lUU]): 

5- _.4 . _J^' (65) 



-43,1 : ei = ^Jo), 62 = 63 = ^^0) 
^3,2 : ei = 5(^ 62 = 63 = 5(^0), 



^3,3 : ei = 62 = 63 = gf^y (66) 

Non-zero commutators for matrices (l65l) and ( 166|) are [62, 63] = 63 and [61, Ca] = 
ea{oi = 2,3) respectively. The algebras of operators (ITTj) corresponding to real- 
izations fl65l) and flUUj) are of the following general forms: 

< j^D - 61, uD - 62, 63 >, < 61, D + 62 + /ite3, 63 > (67) 

and 

< fiD - ii, 62, 63 >, < D + ei + ute2, 62, 63 >, .^g. 

< D + ei + z/te3, 63, 62 >, < ei, ^^62 + G^e3, ^^62 + ^^63 > ^ ' 

corresp ondingly. 

In addition, we have the only four-dimension algebra 

Ai ■ ei = fi-fo)' 62 = 63 = 5(^0), 64 = (69) 
which generates the following algebras of operators ( 1271) : 

<fiD-ei, uD-is, 62, 64 >, < ei, D + 63 + z/te4, 62, 64 >, .^^x 

< + ei + z/te2, 62, 63, 64 > . 
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Finally, it is necessary to take into account the special type of (m+2) - 
dimensional algebras A generated by two-dimension algebras (13H1) . namely, al- 
gebras whose basis elements have the following general form: < fiD + ei + {at + 
\"PxrjXp)e2, Xue2, 62 > where z/, a, p run from 1 to m. The related classifying equa- 
tions generated by all symmetries 0:162, 2:262, ■ ■ ■ ,2:^62 and 62 coincides and we 
have the same number of constrains for /^,/^ as in the case of two-dimension 
algebras A. 

The case a = 1 appears to be much more complicated. The related matrices g 
are of general form ( l54l) and defined up to the general equivalence transformation 
(155|) with arbitrary K"-^. Namely there are seven non-equivalent matrices ( 154|) . 
including g^,g'^ (!56l) . g^ fl60|) and also the following matrices 



9' 













^0 












V 







1 







1 






^0 


1 





(71) 































1 










1 


) 




V 


1 


/ 



In addition, we have fifteen two-dimension algebras of matrices 

^2,1 =< ^(0), 9'^ >, ^2,2 =< gfo), 9(0) >, ^2,3 =< ^fo)' 9^ >, 

^2,7 =< ^^ / >, ^2,8 =< ^fo)' 9^ >' ^2,9 =< ^fo)' 9^ >' (72) 

^2,10 =<^(i), / >> 

A2A =< g\ / >, ^2,5 =< 9\iy 9^ >, ^2,6 =< 9I0) >, 

v42,ll =< X^^'U^l, >, ^2,12 =< -9^\ g^ >, ^2,13 =< ^(2), >, (73) 
^2,14 =< 9\ ^fo) > 

where 

/ \ 

^1°= 1 , ^J2)=^l|,=2 = 

V 1 0; 

Algebras (172|) are Abelian whereas algebras (173ll are characterized by commu- 
tation relations flUTl) . The corresponding algebras ^ are given by equation (lU^ 
and (lUHj) respectively. 

Three-dimension algebras Aa^^ are the algebras A^^^i — As 3 given by relations 
( I65l) . (1661) where matrices (7^ and (7^ are of general form ( |56l) with arbitrary A (i.e.. 
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g^, etc.) and also algebras ^3 4 — A311 given below: 



4 


Cl 


— ^8 ^ — 

— 9 ^ 62 — 

, / ; C2 


^(1) 


! 63 = 


-9(0) 


^3,5 


p-i 


y ) 


Pn 


y(o)' 


^3,6 


' 1 


y(o)' 


y ) 


C3 


y ) 


^3,7 


ei 


= 62 = 




63 = 




^3,8 


ei 


= ^/^ 62 = 




63 = 




^,9 


ei 


= 9fo), 62 = 




, 63 = 


= ^/^ 


^3,10 


: ei 


= ^^ 62 = 




63 = 




^3,11 


: ei 


= 9(0), 62 


= c/ 







Algebras ^3 4 — ^3 g and ^3 7 are isomorphic to A3 1 and ^3 3 respectively. 
The related algebras A are given by equations flU7|) and flUSl) correspondingly. 

Algebra ^3^3 is isomorphic to ^3^3 and so generates algebra f l68l) . 

Algebras ^3^9 and A3 10 are characterized by the following commutation rela- 
tions 



[62,63] = 61 



(74) 



(the remaining commutators are equal to zero); non-zero commutators for basis 
elements of ^3^11 are given below: 



[61,62] =62, [61,63] = 62 + 63. 



(75) 



Using (17^ and fl75|) we come to the following related three-dimension algebras A 
generated by ^39 and ^310: 



< fiD — 262, — 263, 61 >, < 61, D + 26q, + 2z/tei, 6^' >, 

< 6^*61, 6^^*60,, 6q,' >, a, a' = 2,3, a' 7^ a 

and algebras (!77|) generated by A^ n: 

< flD - 26i, 62, 63 >, < 61, 6''*62, 6''*63 > . 



(76) 



(77) 



Finally, four-dimension algebras of matrices ( l55i) are ^4^1 given by equations 
(!69|) and also A^^^^A/^^'^ given below: 



^4,2 
^4,3 
^4,4 
^4,5 



61 = fi-^', 62 = 63 = g\^y 64 = 
61 = ^fo)' ^2 = ^^ 63 = 64 = g^\ 
61 = g^, 62 = 63 = g^, 64 = 

61 = g"^, 62 = 63 = 64 = g^. 



We do not present the related algebras A because all possible non-linearities 
and p will be fixed asking for invariance of equation ([2]) which respect to 
transformations generated by three-dimensional algebras. 
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